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Abstract 

This paper discusses Random Matrix Models which exhibit the 
unusual phenomena of having multiple solutions at the same point 
in phase space. These matrix models have gaps in their spectrum 
or density of eigenvalues. The free energy and certain correlation 
functions of these models show differences for the different solutions. 
Here I present evidence for the presence of multiple solutions both 
analytically and numerically. 

As an example I discuss the double well matrix model with po- 
tential V{M) = -^M^ + f where M is a random N x N matrix 
(the matrix model) as well as the Gaussian Penner model with 
V{M) = ^M'^ — tlnM. First I study what these multiple solutions are 
in the large N limit using the recurrence coefficient of the orthogonal 
polynomials. Second I discuss these solutions at the non-perturbative 
level to bring out some differences between the multiple solutions. I 
also present the two-point density-density correlation functions which 
further characterizes these models in a new universality class. A moti- 
vation for this work is that variants of these models have been conjec- 
tured to be models of certain structural glasses in the high temperature 
phase. 

PACS: 02.70.Ns, 61.20.Lc, 61.43.Fs 
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1 Introduction 



Matrix models are known to be of importance in many diverse areas, for ex- 
ample, in quantum chaos, disordered condensed matter systems, two-dimensional 
quantum gravity, quantum chromodynamics and strings. In the context of 
the one-matrix models, the models that have been studied correspond to an 
eigenvalue distribution on a single-cut in the complex plane where the eigen- 
value density is non-zero [|l|. An obvious generalization is to study a matrix 
model with a more complicated eigenvalue structure. Here I study a class of 
models with a single hermitian matrix but with two cuts for the eigenvalue 
density and point out some unusual features of these models different from 
the single-cut matrix models. One of the important differences observed in 
these models is that they have multiple solutions whose effect appears in 
certain correlation functions. This work suggests the possibility that these 
multiple solutions may be glassy and hence these models may be useful in 
describing certain real disordered glassy systems ref . 0, |^ . 

I study here the matrix model with a double-well potential as well as 
the Gaussian Penner model where similar results are obtained ref. In 
both cases the potential is symmetric about the origin. In addition to the 
usual symmetric solutions I discuss the symmetry breaking solutions. I study 
these solutions in the large N limit and then at the non-perturbative level. 

1 also calculate the two-point density-density correlators using orthogonal 
polynomial methods and methods of steepest descent which strengthens the 
above observations and classifies the models in a new universality class [^, H . 

2 Notations and Conventions 

Let M be a X hermitian matrix. The partition function to be con- 
sidered is Z = / dMe~^'^'^^^^\ The Haar measure is given by dM = 
Uf^^dMiil\i<jdM^fdMj:j^ where Mij = M^f +iMlf; there are A^^ inde- 
pendent variables. V{M) is a polynomial in M: V{M) = giM + {g2/2)M'^ + 
{g3/3)M^ + {g4^/4)M^ + .... TrV{M) and the measure dM are invariant un- 
der the change of variables M — > M' = UMW where U is a unitary matrix. 
We can use this invariance to express Z as an integral over the eigenvalues 
Ai, As, Atv of M. Writing M = UD'U^ where D' = diag{Xi, A2, A^), the 

partition function becomes Z = J dU J^^Uti rfA,A(A)2e-^S-i ^^^'^ where 
= Y\.i<j ~ ^j\ is the Vandermonde determinant arising from the 
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Figure 1: (a). The confining potential (b). The density of eigenvalues 

change of variables. The integration dU is trivial because the integrand is 
independent of U due to the invariance, and gives a constant factor. By 
exponentiating the determinant one arrives at the Dyson Gas or Coulomb 
Gas picture. Thus the partition function is Z = C UiLi dXie~^^'^^ where 
S{X) = NY^iLi ~ 2I]i^j/^|Aj — Ajl, and C is a constant. 

This is a system of particles with coordinates Aj on the real line, con- 
fined by the potential V and repelling each other with a logarithmic repul- 
sion. The spectrum or the density of eigenvalues p{x) = {jj^Y^iLi ^i^ — Aj)) 
is, in the large N limit, just the Wigner semi-circle for a quadratic potential. 
The physical picture is that the eigenvalues try to be at the bottom of the 
well. But it costs energy to sit on top of each other because of logarith- 
mic repulsion, so they spread, p has support on a finite line segment. This 
continues to be true whether the potential is quadratic or a more general 
polynomial and only depends on there being a single well though the shape 
of the Wigner semi-circle is correspondingly modified. For the quadratic po- 
tential V{x) = the density is p(x) = ^^J{x^ — a^) where x E [—a, a], and 
p{x) = elsewhere. The region [—a, a] is said to be the 'cut' where p has 
support. The end of the cut is given by a = J-. See Fig. |l]. 




On changing the potential more drastically by having two wells the den- 
sity can get a support on two disconnected segments ref . 0, § . The simplest 
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example is the potential V{M) = fM^ + fM^ with < 0,g > 0. When 
the wells are sufficiently deep, specifically, when > 2^, the density of 
eigenvalues is given by 
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p{x) = —xJ (x^ — a'^){h'^ — x^), a < X <b or — b < x < —a 



= otherwise (2.1) 

where = — 2y/g] and 6^ = + 2y/g]. The eigenvalues sit in 

symmetric bands centered around each well. Thus p has support on two line 
segments. As \p\ approaches 2y/g from below, a and the two bands 
merge at the origin. For p > —2^, the density is 




otherwise. (2.2) 

The phase diagram and density of eigenvalues for this case is shown in 
Fig. |. 

The simplest way to determine p{z) explicitly is to use the generating 
function F{z) = {jqTr^-^) . F{z) satisfies a Schwinger- Dyson equation 

whose solution is F{z) = ^[V'{z) + a/A] with A{z) = V'{zf - 46(z) and 
6(z) = gz^ + /i + g{j^TrM^) (see ref. 0]). The density p{x) is then deter- 



mined by the formula p{z) = —■^lm^A(z). In what follows I will outline 
the recurrence coefficient method of the orthogonal polynomials to establish 
that there exist multiple solutions which give the same free energy F, F{z) 
and p{z) in the large limit but differ at higher orders. Then I give the 
results for the two-point correlators (also known as the 'smoothed' or 'long 
range' correlators); these tend to different limits as is taken to infinity 
along different sequences (odd or even). This property is similar to that 
suggested in another model of glasses ref. M. 



3 Orthogonal Polynomial Approach 

The partition function Z can be rewritten in terms of the orthogonal poly- 
nomials P„. These are defined as P„ = A" + Cil^A"-^ + . . . + C^h + 
where Cf""^ are constants, and satisfy the orthogonality conditions (P„, Pm) = 
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Figure 2: (a). The double- well potential (b). Density of eigenvalues (c). 
The phase diagram 
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/~^C/Ae-^^WP„(A)P™(A) = hj^m . For example: Po(A) = 1, Pi(A) = 
A + Ci'\ P2(A) = A^ + cf^X + Cq'^ .... Then the partition function in terms of 
the orthogonal polynomials is, ref. [|T0[] , 



TV 



PoiX,] 



Po(Atv) 
Pi{Xn) 



-P/v-i(Ai) P/v_i(A 



N 



(3.3) 



Z can be expressed in terms of the hnS; Z = N\hQhih2 h^-i- For ex- 
ample: in the AT = 2 case, Z = J dXidX2e-^^^^^^-^^^^^\Po{Xi)Pi{X2) - 
Po(A2)Pi(Ai))2 = hohi + hoh = 2\hohi. 



3.1 The Recurrence Coefficients 

The Pn satisfy recurrence relations, ref. jlOl , 

(3.4) 

where P„ and Sn are the recurrence coefficients which depend on the po- 
tential. These recurrence coefficients are central to our analysis, since the 
free energy and all correlation functions can be expressed in terms of them. 
The reason why the P„'s satisfy such a simple recursion equation is that 
/ xPnPn-2G~^^^^^dx = because xPn-2 is a polynomial of degree n — 1 and 
can be expressed as a superposition of P„_i, P„_2, . . . Pq. Thus P„_2 and 
lower order polynomials do not appear on the right hand side of the recur- 
rence relation eq. ( p.4|) . It can be shown that = /i„_iP„. Thus the product 

hohi h^-i = /lo(^o-Ri) (^0-^1-^2) {hoRi...RN) = /i^Pj^~^P^~^ Rn-i, 

hence the free energy F = In Z = In iV! -|- iV In /iq + J2n=i ~ ^) Rn- 

To solve for the recurrence coefficients we have three methods based on 
(i). Integrals, (ii). Recurrence Relations, and (iii). an Effective Potential. 



3.1.1 Integrals 

The recurrence coefficients can be determined by integrals or moments /„ = 
/ dxx'^e~^^^^'^ as described below. The P„'s can be expressed in terms of 
Rn and as P„+i = xP„ - S'„P„ - P„P„-i, with Pq = 1, Pi = (x - Sq) 
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and P2 = {x — Si){x — Sq) — Ri etc. Consider the equation = {Pq, Pi) = 
J dx{x — 5*0 )e~^^*^^'', this results in an integral equation for the coefficient 5*0 
i.e. 'S'o = / dxxe~^^^^^ = As another example the recurrence coefficient 
Ri can be determined as follows. The integral expression for the recurrence 
coefficient Ri is then found from Ri = where = f dxe~^^^^^ = Iq and 
= /c/xe-^^(^)p2^x) = /rfxe-^^(^)(x - Sof = J dxe-^^^'^^x^ - 2xSo + 
Sq) = h — 2/1 5*0 + 5*0. Hence Ri = ^ — (j^Y- Similar expressions for all the 
other recurrence coefficients in terms of integrals can be found. 

3.1.2 Recurrence Relations 

The recurrence coefficients satisfy recurrence relations that follow from the 
identities 

/. J dxe-^^^''^Pn{x)V'{x)Pn{x) = (3.5) 

and 

//. nhn-i = N J dxe-^^''''^Pn{x)V'{x)Pn^i{x) (3.6) 

I follows from the identity = / Pn{x)P'^{x)e~^^^^'^ which holds because 
P^, being a linear combination of P„_i and lower order polynomials, is or- 
thogonal to P„. II follows from the identity / P'^{x)Pn-i{x)e~'^^'^^^ = nhn-i- 

Let us take the following examples: 

(a) V{x) = fx2. 

For this potential the recurrence relation I is = /i(Pn, Pn+i + SnPn + 
RnPn-i) = Snhn, which implies Sn = 0, while the recurrence relation II 
is nhn-i = fiN{Pn, P„ + Sn-iPn-i + -Rn-i-Pn-2) = /iA^/^n which implies 
= /iP„. This determines exactly the recurrence coefficients S'„ = and 
Rn = ^ for this potential. 

(b) V{x) = f a;2 + fxl 
The recurrence relation I is 

= fxSn + g[Rn+i{Sn+i + 2S'„) + P„(2^„ + Sn-l) + S^„] (3.7) 

while II is 

^ 2 2 

— = fiRn + g{Rn-l + Rn + Rn+l) + + S^_i + Sn-lSn- (3.8) 
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Using the initial values for i?o = 0, 5*0 and -Ri we can determine Si, R2 ■ ■ ■ 
using the above recurrence relations for the recurrence coefficients. 

For this potential we now consider the two cases corresponding to V 
having one or two wells. 
(1) The 1-well case: > 0, g > 0. 

It follows that Sn = 0, because all /„ are zero for odd n whenever V{x) is an 
even function of x. Then the recurrence relation I is trivially satisfied. The 
recurrence relation II is 

- = Rnifi + giRn+l+Rn + Rn-l)). (3.9) 

Thus K+i = - ^ - i?„ - We determine i?i = | - (^)' by 

evaluating /2,/i and Jq numerically and evaluate Rn for n > 1 using this 
equation. The result, shown in Fig. ^ suggests that -R„'s lie on a smooth 
curve. This curve is analytically determined as follows: For the large limit 
we set = X and make the ansatz that i?„ is a smooth function of x and 
expand as Rn = R{x) + j^Ri{x) + -^R2{,x) + . . . . Then to leading order in 
N eq. (|3l9| ) implies that x = R{x)[fi + 3gR{x)] which is a quadratic equation 
in R{x) with the solution R{x) = ± a///^ + 12gx]. This fits very well 

with the numerical evaluation of Rn which can be approximated by a smooth 
curve at large N as shown in Fig. |^. 

The generating function for the single well case is given by (see ref. [Q) 

F{z) = dx . ^ (3.10) 
^0 ^z"^ - AR{x) 

On substituting R{x) in the above equation for the generating function one 
gets the same answer as that found by the Schwinger-Dyson equation. Eq. 
p.lO| ) yields the expression eq. ( |2.2| ) for the density of eigenvalues, with a 



single cut. 

(2). The 2- well case, /i < 0, ^ > 0. 
Fig. ^ exhibits a numerical result for Rn in this region of phase space which 
shows that the assumption of a single smooth function describing Rn is no 
longer correct. It suggests the following ansatz for Rn at large N: 

For n < n ( for some n to be determined), we have a 'period-2' structure: 

P ^ / ^" = ^(^) + if^i(^) + • • • for n even , . 

" \ Bn = B{x) + ^Bi{x) + ... for n odd, ^ ' 
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Figure 3: Recurrence coefficients for the single well 
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Figure 4: Recurrence coefficients for the double well 



and for n > n, we have a 'period-1' structure (as for the one- well case): 

Rn = Rix) + j^Riix) + ... (3.12) 



We denote x = ^. Substituting eq. (|3.11|) into eq. (|3.g| ) and equating equal 
powers of 4 we get, for x < x = 



A{x) = - fi^ - igx], 



B{x) = —M + ^fi^-4gx]. (3.13) 



For X > X, substituting eq. ( p.l2| ) into eq. ( |0| ) we get 



R{x) = ^[-fi + ^fi^ + 12gx]. (3.14) 

The above analytical result and numerical figure agree very well. Note that 
X = 1 is the equation for the phase boundary between 1-cut and 2-cut phases 
(see section 2) in the double well region of parameter space. The generating 
function in terms of the recurrence coefficients for the two well case is ( see 
ref. [0 for a derivation ) 

Fiz) = / dx . , if X > 1, 

(z2 - {A{x) + B{x))Y - AAB 







dx — , + / dx — , — if X < 1. 



'0 ^{z'^-{A + B)y-AAB ^x^-4R{x) 

(3.15) 

When X > 1, i.e. in the 2 cut phase, substituting A,B from eq. ( p.l3|) , 

F{z) = ^[fiz + gz^ — gzJ [z"^ — a?) (z"^ — 6^)] which is the same as obtained 



by the Schwinger- Dyson equation. This yields the solution eq. ( |2.1| )for the 
density of eigenvalues. When x < 1, eqs. (|3.13|) and ( |3.14|) leads to the one 
cut solution eq. (p.2|). 



3.1.3 Effective potential method for determining recurrence coef- 
ficients 

Numerically we can evaluate the recurrence coefficient by minimizing an ef- 
fective potential Veff that can be determined from the recurrence relations 
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e.g. for V{x) = ax + + |a;^, the recurrence relations are: 



^ 2 2 

= RnlfJ' + g{Rn+l + Rn + Rn-l + + S^^_i + Sn^lSn)], 

= a + flSn + 9[Rn+liSn+l+2Sn)+Rni2S^n + Sn-l) + Sl].i3.16) 

It is easy to see that if one defines 

°° -n 

n=0 

+ ^(-Rn + 2-Rn-Rn+l) + Cr5'„ + ^S*^ 

+ ^"^n + 9Rn{Sl + + Sn-lSn)}, (3-17) 

then the recurrence relations follow by setting = and = 0. For 
(7 = 0, setting = in eq. ( |3.171 ) and minimizing with respect to Rn 



yields period 1 or period 2 solutions for i?„ as shown in figures 3 and 4. 
However, when a 0, fi < 0, g > 0, (i.e. V{x) has two wells which are 
asymmetric), then numerical minimization of eq. (|3.17|) with respect to Rn 
and Sn yields a solution of the type displayed in Fig. ^ The figure suggests 
that the recursion coefficients possibly become chaotic. However, whether 
they are really chaotic or whether this is only apparently so, requires more 



detailed numerical work |TT|. It would be rather interesting to do this and 
characterize this chaos and also to understand why the recursion coefficients 
are so complicated. At any rate, this suggests that it might be interesting 



to explore solutions of the eq.'s (|3.16D for non-zero but small a, and see 



whether at a = 0, there are alternate solutions ( other than the one with 
Sn = discussed earlier ). This indeed turns out to be the case. Fig. 
^ shows a numerical solution to eq. (|3.16|) with a = 0, obtained by first 
obtaining solutions for progressively decreasing nonzero values of a, and using 
the solution for a particular value of cr as an initial condition to get a solution 
for the next smaller value of a. As seen in Fig. ^, even at a = 0, we can 
get a solution with S'„ 7^ by this procedure. As we show analytically in 
the next sections, there is a large family of multiple solutions that exist at 
(7 = 0. One possible explanation for the apparent chaos in Fig. |^ is that a 
nontrivial mixing between multiple solutions might be occurring. 
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Figure 5: (a). Recurrence coefficients i?„ for the asymmetric 
Recurrence coefficients Sn for the asymmetric double well 
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Figure 6: Graphs of recursion coefficients for a spontaneously broken so- 
lution of the double-well potential with N = 512, /j, = —2,g = 1 (a). Re- 
currence coefficients i?„ and (b). Recurrence coefficients Sn after 100000 
minimization steps from a random start, (c). Orbit in the A-B vs C-D 
plane. 
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4 Multiple Solutions In The Symmetric 
Matrix Model 



The simplest way to understand the existence of multiple solutions which 
appear in these models is to consider the following integral 

for V{x) = ax + + ^x^ with /i < and g > 0, in the limit a ^ 0, V 

has a maximum at x = and minimum at x = = ±y'^. The integral /i 
in eq. ( |4.18|) is zero if we take a = 0, and then take N ^ oo for then the 



integrand is an odd function. Hence 5*0 = for a = 0. However, let us take 
0" 7^ first, evaluate the integral, in the limit N oo and then take the 
limit a 0. The integral can be evaluated using the saddle point method. 
Since for a ^ 0, one of the minimas gives dominant contribution, the result 
is 5*0 = x+ as cr — i> 0^. Thus 5*0 depends upon the order of the limits iV — > oo 
and cr — > 0. 

A more precise way of establishing the presence of multiple solutions in 
these models is by using the recurrence coefficients; the previous sections 
have built up the tools that will be used here. Let us relax the condition 
Sn = 0; we will find that this will be interesting. Consider a period two 
ansatz for both i?„ and Sn- Then 

Rn An, Sn Cn,ioT u cvcu, and 

Rn ^ Bn, Sn-^DnJoT U odd. (4.19) 

As before in the large limit setting An, Bn, Cn and Dn to be equal to 
smooth functions of x = denoted A, B, C , D respectively, one finds that 
eq.'s ( |3.16|) with a = reduce to 



n even x = A[fi + g{2B + A + + + CD)] 

n odd X = B[fx + g{2A + B + + + CD)] 

n even = ^iC + g{B{D + 2C) + A{D + 2C) + C^) 

n odd Q = jiD + g{A{C + 2D) + B{C + 2D) + D'^). (4.20) 

There are four equations and four unknows here but after some work we 
find that there is only three independent equations. The three independent 
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equations are 



C + D = 

A + B + C^ = ^ 
9 

AB = -. (4.21) 

9 

Thus there is a infinite class of solutions labelled by one function of x in 
the large limit. For example, let us consider the two extreme solutions, 
(i). The 'symmetric solution', characterized hj C — D = 0. Then eq. ( [4.21[ ) 



implies C = D = 0, and A{x), B{x) are given by eq. ( p.lll) . This is the same 
solution as discussed in section 3. (ii). The 'maximally asymmetric solution' 
characterized by, A — B = 0. Then eq. ( [4.21| ) implies A = B = R = ^ and 

C = —D = [liii — .f^]^ . The entire infinite class of solutions have the same 
eigenvalue density and free energy in the large limit. This can be seen by 
evaluating the generating function which turns out to be 

l{z^ ~z{C + D)~{A + B- CD)f - AAB 







Eq. ( [4.22|) contains precisely the same three combinations that are fixed 
by the recurrence relation eq. ( ^4.21| ). Therefore independent of which solu- 
tion is chosen we get the same F{z). Since F{z) determines p and the latter 
determines F at large N, this proves that in the limit — ^ oo we have an 
infinite set of solutions of the recurrence relations with the same eigenvalue 
density and free energy. This demonstrates the presence of multiple solutions 
from the recurrence coefficient point of view. 

5 Non-perturbative Solution 

The multiple solutions found above show differences at higher powers of j^, 
for example, in the double-scaling limit ref. 1^, 1^, 14, 15, 16, 17, 18 1 which 



I describe in this section. Let us begin by taking the symmetric solution and 
proceed by expanding the even, odd recurrence coefficients A, B as 

A„ = ao + e{f^{t) + foit)) + e^{r,{t) + ro{t)) n even 

Bn = ao + e(/e(t)-/o(t)) + e'(re(t)-ro(t)) n odd (5.23) 
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where x = ^ 



1 — eH, ao -- 
powers of e, we get fe = 



^2^, and e = N Then upon equating equal 
0, /o = fit), Te = where the susceptibihty 



is given by x 



. f satisfies f" — \P + \ft = the Painleve II equation. 



We then take the asymmetric solution with A 
making the following expansion 



Cn 



eg{t) + e\.. 
-eg{t)+e\.. 



B ^0, C - D ^0. On 



(5.24) 



we get the following coupled equation ref. [jT^ [T^, |T7|, [TB| /" — \f{g^ + 
P) + ^ft = and g" — \g{g'^ + P) + \gt = 0. If we make the substitution 
/ = rcosO, g = rsin^ where r and 6 are functions of t, then r satisfies a 
modified Painleve equation 



1,1 /3 

-r H — tr - 

4 2 r3 



0, r 6 = I = const. 



(5.25) 



4 



The case Z = is the Painleve II equation. The susceptibility x 

and X ~ T ~ (^~4^)^~^ + • • • which is very different from the symmetric 
case. Thus the multiple solutions show differences at higher order in j^. 



6 Distinguishing multiple solutions: Correla- 
tors, Odd and Even N 

Multiple solutions are distinguished by certain correlators. Consider the cor- 
relator (TrM TrM)e = (TrM TrM)-{TrMf where (O) = | / rfMe-^^^^(*^)0. 
It can be shown |^ that {TrM TrM)^ = Rn- Since Rjy depends on the 
choice of the solution, this demonstrates that the multiple solutions under 
discussion give rise to different correlation functions in general. In particular 
note that -Rat = if N is even and Rn = if N is odd. For the 'sym- 
metric solution' A ^ B, hence this correlator changes by 0(1) as goes 
from odd to even. However for the maximally asymmetric solution A = B, 
hence this correlator remains essentially the same as is changed from odd 
to even. 

Another example of a correlator which distinguishes between the two 
solutions is 

(TrM TrM TrM)^ = Rn{Sn^i - Sn). (6.26) 
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For the symmetric solution S'„ = for all n, hence this correlator vanishes. 
For the asymmetric solution S'„ is period two and the correlator changes sign 
while going from odd to even A^. 

The orthogonal polynomial P„ can be thought of as a 'wavef unction' of 



a 'state' \n) in the 'coordinate basis', i.e. (A|n) = J^-^ e 2^y^\ where |A) 
are eigenstates of the operator M. The operator M is defined in the basis of 



states \n) as the matrix 



5*0 

/R'l 





yRi 

Si V R2 

/ R2 5*2 

which follows from the recurrence relation satisfied by the orthogonal poly- 
nomials eq. 



see ref. M 



(6.27) 



for more details. It can be shown that the 
density of eigenvalues of this matrix in the large N limit is precisely p(A). 
It is therefore of interest to calculate the eigenvalues of this matrix. Since 
the matrix elements are given by the recurrence coefficients we can deter- 
mine the eigenvalues numerically (see section ( |3.1| )) from a given solution of 
the recurrence coefficients. The figures show the location of the eigenvalues 
for the symmetric solution and the maximally asymmetric solution (see Fig. 

For N even, half the eigenvalues are in one well and the other half in the 
other well. While this is true for both solutions the detailed positions are not 
the same. When is odd, for the maximally asymmetric solution one extra 
eigenvalue is located in one of the wells (the well it is located in depends 
on the sign of 5*0). For the symmetric solution this extra eigenvalue is in 
the center, thus preserving the symmetry between both wells. This seems to 
suggest that the bulk effect of the two solutions is the same but they differ 
by one eigenvalue effects. 

Further evidence for the effects of multiple solutions comes from corre- 
lation functions of the density operator p{x) = j^TrS{x — M). The explicit 
expression for the smoothed two-point correlator for a general 2-cut solution, 



obtained by the method of steepest descent [|^, ^ is 
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m 




Figure 7: Eigenvalue distribution for even and odd N: (a) N=24, symmetric 
solution; (b) N=24, asymmetric solution; (c) N=25, symmetric solution; (d) 
N=25, asymmetric solution. 
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Here a{z) = Y[i=i{z — s = ai + 02 + as + 04, (the support of the eigen- 
values consists of the two segments [01,02] and [03,04]), eA = +1 for AG 
[03,04], eA = — 1 for A G [01,02] and (3 = 1,2,4 depending on whether M 
the matrix is real orthogonal, hermitian or self-dual quartonian. C is an 
undetermined constant in this method. It turns out that the same correlator 
can be calculated using other methods, which yield different values of C. For 
example in [^ we obtained an asymptotic form of Pn{x) for the symmetric 
double well potentials for n close to A^. Using this form of P„ (which corre- 
sponds to the 'symmetric solution'), we found the smoothed < pp >c given 
by eq. ( |6.28| ), with C = (—1)^, (the origin of this nontrivial N dependence, 
also observed in |2^, is explained in [^). On the other hand in ref. [24 



this correlator was calculated by the loop equation method starting from an 
asymmetric V and taking the limit of a symmetric potential where C was 
found to be — |{(o^ + 6^) — (a + &)^[^^]} where E{k) and K{k) are complete 

elliptic integrals of the first and second kind and k = We conjecture 

that different values of C corresond to different solutions of the recurrence 
coefficients. 

In this section, we have shown the multiple solutions discussed in section 
4 give rise to different correlation functions. An interesting feature of the 
solutions is that there is a nontrivial N dependence which survives in the 
large N limit. In particular, for the symmetric solution we have shown a 
difference between odd and even in the thermodynamic limit. This is 
reminiscent of the behaviour in a model for glasses 



7 The connection of random matrix models 
with a simple model of structural glasses 

This section is concerned with building up a connection with the study of 
structural glasses in the high temperature phase, starting from the work of 
ref. [0, ^. The hamiltonian H = ^J^aftbi^'^-^'^Y corresponds to n particles 
moving in an A^-dimensional space, in the limit n, N 00. The coordinates 

of the particles are = (s" s%) and o = 1 n. In the 'spherical' case 

the particles are constrained to move on the surface \s"'\'^ = N, for all o. In 
the 'Ising' case they occupy only the vertices of a hypercube = ±1. 

The corresponding partition function is Z = e^^ Tr[s}Gxp[—-i^tr{S'' SS"^ S)] 
where S is the N x n rectangular matrix of elements s", and Tr^s} runs over 
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either the spherical or the Ising measure. Here n = aN with a > 1. Using 
the large N equivalence with the global constraint one gets 

- (3E[x]j (7.29) 

The Xi\fN are the "diagonal" values of S in its canonical form and E[x\ = 
NY^f=i y{p^i)~^ ^i^^j ^iT'lxf—x'jl where x = {xi, X2, xn) and V{x) is given 
by 

This potential has the form of the generalized Penner model. I show briefly 
(see ref. 0) that a very closely related model that can be exactly solved, 
the Gaussian Penner model, has 2 cuts and multiple solutions. So the above 
model which is a variant of the Penner model with two cuts, is also conjec- 
tured to have multiple solutions. 

The potential for a general Penner model is V{M) = Vq{M) — tlnM, 
where Vq is a polynomial. If Vo(M) = |;uM^ the model is the Gaussian 
Penner model where we can re- write InM = ^ InM^. This has Z2 symmetry 
and the potential is a double well with eigenvalues distributed in disconnected 
segments. The partition function for the Gaussian Penner model in terms of 



its eigenvalues Xi is [25 



Z = jY[dxiexp[-E{x)] (7.31) 
where E[x] = N Eti Vi^i) - 2 Ei^j In \xi - Xj \ 



V{x) = N (^x^ -t\n\x\) . (7.32) 



This potential is similar to that of eq. ( |7.30|) , in that it has a double- well po- 
tential involving ln\x\. I will show explicitly that this has multiple solutions. 
We consider the situation for t > 0. 

The recurrence relations for a general Penner model reduce to 

^ = ^n{n-l\V^iM)\n) -t^n{n-l\M-^\\n) 
= {n\Vo{M)\\n) -t{n\M"'^\n) (7.33) 
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Denoting Wn = y/R^{n - l\V^{M)\\n) and r„ = {n\V^{M)\\n), for the Gaus- 
sian Penner model Wn = fiRn and Yn = fiSn- We can consider as in section 
4, a period-2 ansatz for Rn and Sn which leads to four equations but again 
only three independent equations: 



C + D = 
A+B-CD = ^ 

AB = (7.34) 

Thus there is an infinite class of solutions labelled by one function of x 
in the large N limit. For the 'symmetric solution': C = D = and A 



X 



B = while the maximally asymmetric solution: A = B = j^^Jx{x + t) 



and = ^[(2x + t) — 2^x{x + 1)]. Once again eq. (|7.34|) fixes the same 
combinations which appear in the generating function of the Gaussian Penner 
model. Thus in the large N limit the eigenvalue density and free energy are 
identical for the full infinite class of solutions satisfying eq. (|7.34| ). 



For the symmetric solution (n|M~^|n) = by Z2 symmetry and Sn = 0, 
thus eq. ( [7.33|) yields Rn = for n even and i?„ = ^ + ^ for n odd. 
This is an exact solution, hence the exact free energy may be found to be 
r = fclog[(2fc + + l)(2fc + /i - 1)], where t = -1 + f . Expanding 

in powers of /i we get 

r = ^/^^log/x+ ^log/i.... (7.35) 

Note that the coefficient of the second term log /i is Xi = 3^ which corresponds 
to the first subleading correction in the expansion. For the maximally 

asymmetric solution in the double scaled limit Rn ~ r(i/2(jv-n+^+i/2)j ' ^^"^ 
the free energy is P = eS"^ ^ log[(2fc + /i + 1/2) (2A; + /i - 1/2)] + . . . . On 
expanding in powers of /x the free energy is 

1 5 
r = |/x^log/i - — log/i. . . (7.36) 

The coefficient of the second term is 4^. Here in the Gaussian Penner model 
we see that though the symmetric and maximally asymmetric solutions give 



21 



the same answer in the large hmit, the free energies are very different at 
higher order. This estabhshes that in the Gaussian Penner model multiple 
solutions are present which give the same free energy in the large limit but 
differ at higher orders. As in section 7, certain correlation functions will be 
different for these solutions since they depend upon the recursion coefficents. 

The potential V{x) for the glass model discussed in the high temperature 
region is qualitatively similar to the V{x) for the Gaussian Penner model 
in that both have double wells. This feature is the same for the model 
discussed in earlier sections. Since these later models which have two cuts 
in their eigenvalue density have been shown to have multiple solutions, it is 
conjectured here that the glass model described above in the high tempera- 
ture phase should also have multiple solutions. It would be useful to show 
this explicitly in the future. 



More recent work suggests that the number of multiple solutions 
grows exponentially with A^. This gives further support for glassy behavior 
in these gapped random matrix models. 



8 Conclusions 

To conclude, ample evidence has been provided for the existence of multiple 
solutions in random matrix models with gaps. First a simple motivation 
for the presence of multiple solutions is given and then made more precise 
using the recurrence coefficients of the orthogonal polynomials of the system. 
Further, numerical evidence is given for the existence of multiple solutions 
in this context. In the large A^ limit the free energy, generating function and 
density are the same. Differences between the multiple solutions are seen in 
the free energy at higher orders in -j^ as well as in correlation functions. 

Connections with the high temperature phase of structural glasses have 
been made to matrix models ^ . These are a variant of the Penner model 
with two-cuts. A simpler model, the Gaussian Penner model, with discon- 
nected segments is shown to have these unusual multiple solutions as well. 
Hence the matrix models with gaps with connections to glass models are 
likely to have this unusual property. The ruggedness of the landscape needs 
to be studied. It would be nice to be able to cast these models in the replica 
framework but this remains a difficult task at this point, as the Hubbard- 
Stratanovich transformations which are technically needed for the Gaussian 



random matrix models ref. ||2^ are not available here for the simple M 
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and the Gaussian Penner model or any other gapped random matrix models. 
This is a future goal in this problem. 
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